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Abstract 



We consider the ^(v^i + ^t) ^lodel on a d- dimensional Euclidean space, where all but one of 



Qh, the coordinates are unbounded. Translation mvariance along the bounded coordinate, z, which 

D , 

(— I 

' lies in the interval [0,L], is broken because of the boundary conditions (BC's) chosen for the 
^ ! hyperplanes z = and z = L. Two different possibilities for these BC's boundary conditions 
are considered: DD and A^A^, where D denotes Dirichlet and N Newmann, respectively. The 
renormalization procedure up to one-loop order is applied, obtaining two main results. The first 
is the fact that the renormalization program requires the introduction of counterterms which are 
surface interactions. The second one is that the tadpole graphs for DD and A^A^ have the same z 
dependent part in modulus but with opposite signs. We investigate the relevance of this fact to 
the elimination of surface divergences. 



1 Introduction 

In this paper we consider an interacting quantum field theory model in the presence of boundaries. 
We shall assume that the system is finite along one dimension z G [0,L], and infinitely extended 
along the remaining (rf — 1) directions. 

The presence of geometric restrictions on the domain of one of the coordinates of the system, 
demands the introduction of classical boundary conditions, to be satisfied by the fields on the 
two hypersurfaces at z = and z = L. If we restrict ourselves to a real scalar field, Hermiticity 
of the Hamiltonian leads us to five different (inequivalent) choices for the BC's, namely: DD, 
NN, DN, periodic and anti-periodic. The last two choices are usual in the finite-temperature 
literature, and shall not be dealt with here, since they don't break translation invariance, which 
is the phenomenon we are concerned with. 

Physical systems will be, in general, finite along several directions. For the sake of simplicity 
we will consider a d-dimensional layered geometry. Although the highly idealized case of planar 
boundaries misses a whole series of features that are present in the general, curved boundary 
case, for more general shapes the multiple reflection method can be used to find the correlation 
functions of the model [0 . 

Most of the papers in the literature deal with periodic or anti-periodic boundary conditions, 
where translation symmetry is maintained, and surface effects avoided. Moreover, in quantum 
systems where translation symmetry is broken, the renormalization procedure is more involved 
than for translation invariant systems, either bounded or unbounded. 
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The diagrammatic expansion and the renormahzation program for an unbounded system is 
conveniently performed in momentum space. On the other hand, when DD or A^A^ BC's are 
implemented, one may still work with Green's functions at fixed {d—1) dimensional momenta, since 
there is translation symmetry along those dimensions. As discussed by many authors, associated 
with the breaking of translational invariance, a new feature emerges: the existence of one particle 
reducible primitively divergent diagrams. For example, the bare two-point function Gq{x, x') of the 
scalar model with zero, one or two points on the surfaces have different renormahzation constants, 
respectively. 

In this paper we will consider a scalar theory subject to two different classical BC's: DD and 
A^A^. Besides the lack of translational invariance, we shall face the problem of surface divergences. 
One way to avoid them is to smooth out the plates surface. But in this case an ambiguity appears, 
since loop-graphs will depend on an ad-hoc model assumption, namely, the particular features of 
the smooth walls. Consequently, we prefer to maintain the hard walls assumption. In the context 
of the Casimir energy of minimally coupled scalar fields, many authors used soft, hard and semi- 
hard BC's 0. Different questions sometimes require more complicated BC's, like the quantum 
mechanical treatment of the boundary conditions presented by Ford and Svaiter 0, a device 
implemented to solve a long standing paradox concerning the renormalized energy of minimally 
and conformally coupled scalar fields. 

Besides the above mentioned effects due to the existence of surfaces and the breaking of trans- 
lation symmetry, we do also have, of course, finite size effects, which are of a different nature. 



Various investigations have been made on this subject, mostly from the quantum statistical me- 
chanics point of view, and we present a short review here. Pathria and co-workers studied the 
Bose-Einstein condensation of an ideal relativistic Bose-gas confined to a rectangular box of sizes 
Li, L2 and L3 with periodic boundary conditions on all the walls |^. A systematic study of finite 
size systems and phase transitions was developed by Brezin and Zinn- Justin [^. These authors 
studied the 0{N) model in two different geometries: the periodic cube and the cylinder along one 
dimension (the time) and finite and periodic in the {d — 1) remaining dimensions. Nemirovsky 
and Freed considered the same model but in a 'complementary' situation regarding the boundary 
conditions, namely, {d — 1) dimensions are unbounded, and periodicity along the only finite di- 
mension, and p. Afterwards, Singh and Pathria studied the 0{N)d model confined to geometries 
with periodic boundary conditions in all directions [0. The same model in the presence of one 
mirror localized ai z = was also studied by many authors 0. 

Another ingredient, important in finite size systems at criticality, is the concept of finite size 
scaling. Let us consider a finite system of linear size L and suppose a thermodynamic quantity 
Phit) (where t is the reduced temperature) becomes singular a.s t -^ 0. Defining -^rJ- = g{L,t), 
where P{t) is the bulk value of Piit), the statement of finite size scaling is that g{L, t) = f{L/^{t)), 
where ^(t) is the correlation length. In other words, finite size scaling predicts that, for large L, 
the dependence on L of the singular contributions to thermodynamics functions scales with the 
correlation length, and is described by universal scaling functions. It may be pointed out that 
finite size systems must be classified in two distinct groups with respect to finite size scaling. 



depending on whether the B.C.'s break translation invariance or not. For finite size systems where 
translation invariance is maintained (for example, a periodic cube, or a cylinder infinite along one 
dimension, and finite and periodic in the {d — 1) other dimensions), finite size scaling is easily 
understood. 

The renormalization group equations are insensitive to finite size effects (since the renormaliza- 
tion is related to short distance singularities) and must, accordingly, be maintained in such finite 
size geometries. However, the solutions to those equations must be different from those for the 
unconstrained systems, because correlation functions can depend on the additional dimensional 
quantities (the lengths of the compactified dimensions) and finite size scaling is present. 

For the cases of DD, NN or DN B.C.'s, the situation is quite different, since it is much more 
difficult to decide if a given interaction is relevant, irrelevant or marginal, the reason being that 
the propagator of the critical theory satisfies B.C.'s which can interfere with the power counting. 

Finite size effects have also been extensively studied in the quantum field theory context 
during the last twenty years. In fiat spacetime with one compactified dimension, the mass can 
depend upon the periodicity length 0. This phenomenon is of particular interest in theories with 
broken symmetry, as it allows topological effects to play a role in the restoration of symmetry. 
An equivalent mechanism is at work when we assume that the fields are in thermal equilibrium 
with a reservoir at temperature (3^^ |Tn|. Finite size effects in quantum field theory with periodic 



boundary conditions in the spatial section at finite temperature was analyzed by many authors. 



See, for example, [11|, and references therein. 



In the context of non-Abelian gauge field theories at zero temperature, cavity QCD was studied 



by many authors |T2|. Hanson and Jaffe and also Hanson, Jonhson and Peterson dealt with 
quantum fields in bounded domains with broken translation invariance. Regarding Abelian gauge 
theories, quantum electrodynamics in the presence of conducting plates has also been the focus 



of research |T3|. Based on the fact that the B.C.'s for the electron field could lead to additional 
contributions to the Casimir force, Bordag et al. and also Robaschik et al. adopted the following 
model. A photon field obeys classical B.C.'s on perfectly conducting plates, while the B.C's for 
the fermion field are free. These authors assumed that the electromagnetic field also exists in the 
region outside the plates (two simple connected domains). Chodos and Thorn investigated the 
self-energy of fermions used different B.C's (the slab-bag), where the fermionic field is confined 
between two parallel plates and the photon field is unconfined. In fiat spacetime, for systems where 



some dimensions are compactified but translational invariance is maintained, Toms WM and also 



Birrel and Ford |15| have shown that all the counterterms are independent on the compactified 



spatial size. A more general discussion has been given by Banach ||I6[. This author proved that a 
topological identification (periodic or anti-periodic B.C.'s) does not introduce new counterterms 
into the theory. As stressed by many authors, were this not the case there would be a catastrophe 
in the renormalizability of the model. 

For translation invariant systems, because of Poincare invariance, one should expect that over- 



lapping divergences will not obstruct the implementation of the renormalization program |T^. In 



systems where Poincare invariance does not hold, these proofs do not apply, and one must show 



that it is still possible to iinplement such program. A technical difficulty is also met here, since the 
presence of geometric restrictions makes Feynman diagrams harder to compute than is ordinary 
quantum field theory in unbounded systems. An crucial work on this subject has been presented 
by Symanzik [TE\. Of particular importance are also the papers by Nemirovsky and Freed, and 



Krech and Dietrich [p!9|j . 

In this work we shall consider an anisotropic scalar model, in a d-dimensional Euclidean space, 
where the first {d — 1) coordinates are unbounded and the last one lies in the interval [0, L]. We 
analyze two different translation symmetry breaking B.C.'s: DD and A^A^ on the plates. We first 
present a rederivation of the fact that to renormalize the theory one has to introduce counterterms 
as surface interactions. We also show that the tadpole graphs for DD and for A^A^ B.C.'s have 
the same modulus for their 2;-dependent parts, but their signs are opposite. We study the possible 
use of this property to get rid of the surface divergences. 

The organization of the paper is as follows: In section II we present the general formalism. In 
section III we discuss the slab configurations, dealing with the two-point and four-point functions, 
both for DD and A^A^ b.c. In section IV we analyze the divergences of the translational invariant 
part of the tadpoles. Section V deals with the analysis of the ultraviolet and infrared divergences 
of the z-dependent part of the tadpoles. Finally, section VI contains our conclusions. Throughout 
this paper we use h = c = 1. 



2 General formalism and the scalar anisotropic model 



Let us consider Z[J]^ the generating functional of complete Green's functions for a scalar field 
in a (i-dimensional Euclidean space 

Z[J] = /r'v^e~'^[^i+/-^(^)^(^) (1) 



where 

S[^] = jd''x[C{^{x),d^{x))], (2) 

Dip is the appropriate measure, and S[ip\ is the classical action associated with the scalar field. 
The quantity Z[J\ can be regarded as a functional integral representation for the imaginary time 
evolution operator ((/92I U{t2, ti) |v?i) with the boundary conditions: v^(ti, x) = ipi{x) and ip(t2, x) = 
(f2{x). The quantity Z[J] gives the transition amplitude from the initial state \ipi) to the final state 
\ip2) in the presence of some scalar source J{x), of compact support. Regarding the Lagrangian 
density C, we shall assume it to be 

C{^) = lid^r + lmV + l\v'- (3) 

The n-point correlation functions are given by the expectation value with respect to the weight 
^-sif) defined as 

G^-\x,,X2,..,x^)=<^{x,).M-^)> = ^^JJl^^^^^JJ^^h-0 



V^^{xi)...^{xn)e~'^''^ . (4) 

As usual, W{J), the generating functional for connected correlation functions of the elementary 
fields shall be given by W{J) = lnZ{J). Thus 

G^^\xi,X2,..,Xn) = , . , . , r|j=0 =<^{Xi)...^{Xn) >c ■ (5) 

J [XijO J {X2) ■ -0 J [Xn) 

Finally, r((/9o), the generating functional of connected one-particle irreducible correlation functions 
is introduced by performing a Legendre transformation on Vr[J], 



r(^o) = -WiJ) + J (fx ip{x)J{x) (6) 

and 

Y(n)( ^ N S^^M I .^x 

dipo{Xi)d(po{X2)..dipo{Xn) 

where 

^»<^> = jwr '"' 

If A = the partition function Z{J) can be calculated exactly i.e. 

Zo( J) = ea;p Q | rf'^a^rf'^y J(a;)Z}(a; - y, m')J{y)^ , (9) 

where 

{-A.^ + m^)D{x-y,m^) = 5'^{x-y). (10) 



For A 7^ it is not possible to find exactly Z{J) and perturbation theory is mandatory. This 
expansion stems from the formal identity: 



Z{J) = exp 



A\J \6J{x) 



ZoiJ). 



:iii 



From now on, we shall consider a generalization of the previous case, namely, the anisotropic 
Landau-Ginzburg model for a N = 2 component order parameter with a Lagrange density C = 

^0 + Qnt, where 

1 „ 1 _ 1 „ „ 1 _ _ 

(12) 



1111 
Alv'i, V^2) = 2 (^V^i)^ + 2^9<^^f + 2"^V? + 2"^V2 



and 



^int = ■^iVl+Vt)- 



(13) 



To generate the n-point functions we have to introduce two scalar sources Ji{x) and J2{x) coupled 
linearly with the fields ^pi{x) and ip2{x) respectively. Integrating out the fields, we obtain Z{Ji) 
and Z{J2) and the total partition function of the model factorizes: Z{Ji^ J2) = Z{Ji)Z{J2), where 



Z{J: 



1,2 j 



Ni, 



-exp 
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'A\I'^'''UJ.2{x) 



Zoi-J: 



1,2 j 



(14) 



and 



^0(^1,2) = exp {^j d''xdS.h,2{x)Gfl{x - y,m^)J,,2{y)) , • (15) 

In the above equation, Gi {x — y,m'^) and G2 (x — y^rn^) are the free propagators, solutions of 
the inhomogeneous equations {i = 1,2). 



.2\M^) 



{-A^ + m')Gr{x-y,m') =5''{x-y). 



(16) 



The partition function applies to arbitrary geometries, and classical B.C.'s must be implemented 
on the Green's functions. As discussed before, we will assume that the system is confined between 
two parallel plates localized at z = and z = L, and use the Cartesian coordinates x^ = (r, z) 
where r is a (d — 1) dimensional vector perpendicular to the z direction. A question that arises 
in such a model is related to the renormalization conditions. It is well known that for fields inter- 
acting with a thermal bath defined in manifolds where the spacelike sections are non-compact the 
mass and coupling constant counterterms are temperature independent. Using dimensional regu- 



larization ||20[ it was proved that for fields defined on manifolds where the spacelike sections are 
non-compact, or compact in at least one dimension, but with the other dimensions noncompact- 
ified, the mass and coupling constant counterterms are size and temperature independent at the 
two-loop level. In a perturbative scheme, the renormalized theory is fixed by the renormalization 
conditions for the superficially divergent vertex functions (the one particle irreducible parts of the 
connected Green's functions). In other words, in the conventional renormalizable (translational 
invariant) theory the ultraviolet divergences can be absorbed by counterterms related to the field, 
mass and coupling constant. A question of fundamental importance is how the renormalization 
program can be implemented in systems where translational invariance is broken. The purpose of 
the next section is to analyze this question for the case of the anisotropic model at the one-loop 
approximation. 
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3 Finite size effects and classical boundary conditions 



For the cubic anisotropic model, we define the boundary conditions over the plates for the 
fields fi(x) and <f2{x). For the (fi{x) field we assume Dirichlet-Dirichlet boundary conditions i.e: 

(pi{f, z)\^=o = (pi{f, z)\^=L = 0, (17) 

and for the ip2{x) we will assume Newmann-Newmann boundary conditions, i.e.: 

d d 

-Q-^2{r,z%=Q = -^'f2{r,z)\^=L = 0. (18) 

It is well known that '^He films close to the A transition satisfies DD b.c. Another well known 
example of such kind of boundary conditions is the electromagnetic field. It was shown that 
for an electromagnetic field confined in a perfectly conducting cavity, it is possible to treat the 
electric and magnetic modes separately, where each one satisfies Dirichlet and Newmann B.C.'s, 
respectively |^ . Going back to our discussion, since the translational invariance is not preserved. 



let us use a Fourier expansion of the fields in the following form: 

(ztt) 2 "' „ 

where the Un{z) are the normalized eigenfunctions of the operator — ^ satisfying the completeness 
and orthonormality relations, i.e.. 



E 



u^{z)u:iz')=Siz-z'), (20) 
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and finally 



L 

dzUn{z)ul,{z) = 6n,n', (21) 



— M„(z) = klun{z), (22) 



where fc„ = ^, n = 1, 2.. for DD b.c and n = 0,l, 2.. for A^A^ b.c. 

It should be noted that this kind of expansion in an orthonormal set corresponding to the 
eigenfunctions of the Hermitian operator — ^ defined on a finite interval could be quite straight- 
forwardly generalized to different anisotropic models. For example, we might consider a scalar 
field defined on Euclidean space with all the coordinates unbounded, but with a mass having an 
anisotropy along the z coordinate. Namely, the Lagrangian density could be 

^ = lid^r + li^" + /^' W)V^' + ^^' ' (23) 

which is, of course, non invariant under translations in z, except for the trivial case of a constant 
fj,{z). If now we assume the Un{zys to denote the normalized eigenfunctions of the Hermitian 
operator 

with 

hUn{z) = \lUn{z) , (25) 



the expansion (|T^) still holds. We assume /x to be a non negative function, so that h is definite 
positive. 
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We see that the case of Dirichlet B.C.'s could be obtained starting from an unbounded z 
coordinate, and using the anisotropic mass /i(z) = /i(|;)", with ra — > cxd, and /i a (positive) 



constant. For a study of Dirac fermions with a space dependent mass see for example ref. p2 |. 

Coming back to the case of DD and A^A^ boundary conditions, the eigenfunctions are, respec- 
tively, 

/2 TITC Z 

-sin(-^) n=l,2.. (26) 



and 



/2 TITT Z 

_cos(-^) n = l,2,--- . (27) 



For A^A^ B.C.'s, we also have the zero mode Uo{z) = -^. The free propagator can be expressed in 
the following form: 



G^o\r,z,z') = -1^ f d''-'pY.e'^^-''ur,iz)u:iz')GoAP), (28) 

where is not difficult to show that G'o,n(p) is given by 

Go,n{P) = ip' + kl + m'r\ (29) 

For the anisotropic mass case, we would have instead: 

GoAP) = ip' + >^l + m'r'. (30) 

As we discussed before, for translational invariant systems we have Gq (x, x') = Gq (x — x') 
and from coordinate space Feynman rules we can go to momentum space representation, which is 
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the more convenient framework to analyze the divergences of the theory. Translation invariance is 
reflected in momentum conservation conditions. Since our system possesses translation invariance 
along the direction parallel to the plates, the parallel momentum, p, is conserved. In this case 
a convenient representation is a mixed {p, z) space. The Feynman rules for different boundary 
conditions was derived in many references and we will not repeat the rules here. For a careful 



study of Feynman rules in such systems see for example Ref. [^. Let us study the one-loop 
correction to the bare two-point function Gq (x, x'), both for the DD and A^A^ cases. Using the 
Feynman rules (see fig.(l)) we have: 



Go (A,ri,2;i,r3,23) = - j (f V2 / 6^22^0 (n, ^1; ^^2, ^2)^0 (f2, ^2; ^2, ^2)^0 (^2, 2^2; ^5, 2:3). 



(31) 



or 



Gj'HA,ri,Zi,f3,Z3) = \j d'^-^r^jl^ dz2G'i\n-f2;zi,Z2)Gf\Q,Z2)Gt\r2-f:,-Z2,z^). (32) 




{n,zi) {r2,Z2) {r3,Z3) 

Figure 1: The two point function for (fi{x) 
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Each of these expressions can, for the case of DD B.C.'s, be expanded as 



and finally 


- ^ ^ ) - 2 1 V ^inf ^^"^ ) ^inf^^"^ ) 1 ^'-\ ""^"'"^ 


(33) 
(34) 

(35) 








^" ^°'^^^-L(2.)<^-^^^/^" L i^ ^(p-'2 + (n;.)2+^2)- 



Although the functions Gq (ri — r2, Zi, ^2) and Gq '{r2 — "^3, -22, -^s) are singular at fi = r2, -^i = -22 
and r2 = "Hs, -22 = -^3, the singularities are integrable (for points outside the plates), consequently 
only the tadpole is divergent and needs a regularization and renormalization procedure. A straight- 
forward calculation yields the order A correction to the bare two-point function in the one-loop 
approximation : 

^(2),, ^ ^ ^ 2 1 /-i ^ . rnizi . 7171X2. . ,n'nz2. . .n'nzs 

Gy{X,ri-r3,Zi,Z5) = JIWY^ J ^ sin(— — ) sin(— — ) sin(— — ) sin(— — ) 

V / n,n'=l 

where, since we will use dimensional regularization, we introduce a dimensional parameter /i, and 
define g = A/x'^"'^. The expression for the tadpole TooiL, m, d, z) is then: 

Too(L,,n,,,^) ^ |^|.„^(H£)/,-V^^,-^A__, ,37) 

The tadpole graph in the case of N-N B.C.'s can be also easily found, and it is given by 

Q \ [' 1 

Tnn{l, m, d, z) = j-pTTTi / ^'^'^^TrT, — ;t 

L (27r)« ^ J (A;2 _^ ^2) 
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^t^^^^'Jll, I ,^-^ I ,3, 



2g 1 ^ _^,2,n7vz 



Note that both Tdd{L, m, d, z) and Tnn{L, m, d, z) diverge in their continuum momenta integrals 
and also in the n summation. In the next section we will analyze the ultraviolet behaviour of 
the bare two-point functions i.e TD£i(L,m, (i, z) and TNN{L,m,d,z) . Before dealing with the 
renormalization of the one-loop two point function, let us, by the sake of completeness, discuss 
the bare four-point function. The expression for the bare four-point function is given below, 
and the same analysis of the divergences can be done. In this paper we will not implement the 
renormalization program of the four point-function (which follows using the same procedure used in 
the two-point function). Our object of interest is the two-point function, since it is the fundamental 
quantity that measures the vacuum activity. Using the Feynman rules, Gq {\,Xi,X2,Xs,X4), the 
order A^ correction to the bare four-point function, is given by 

G^o\X, fi, zi,f2, Z2, fs, Z5, fe, zq) = -j d'^~^r^ j d'^'^r^ j dz^ j dz^ G^Q\fi - fg, zi, z^) 
Gf\f2 - fg, Z2, z^)[G^i^\r^ - n, Z3, Zi)fG'i}\n - ^5, Z4, z^)Gf\n - rg, Z4, Zq). (39) 

Again, all Go's are singular at the same points, but the singularities are integrable, except for 
Gq (r*3, Z3, r4, Z4), consequently, to renormalize the bare four-point function we have the regularize 
the IPI four-point function r*^^)(A,r3 — r4, Z3, 2:4) = [Gq (A, rs — r4, Z3, ^4)]^, which, for DD B.C.'s, 
is given by 

Go (r3-r4,Z3,^4)--^^^^^l.sm(-^)sm(-^)yrf ^(p2 + (-)2 + ^2) • (^0) 
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{r2,Z2) 



{n,zi) 




{r5,Z5) 



[re,ze) 



Figure 2: The four point function for ipi{x) 

A convenient way to express 6*0(^3 — "^4, -23, -24) is the following. Let us define p = r-^ — f^, and also 
Z3 — Z4 = u, Z3 + Z4 = V then it is possible to write 



.(2) 



(2), 



.(2), 



G'o\fs - u, zs, z,) = G'o\p, u) + G',\p, v) 



(41) 



where 



and also 



G'i\p.u) 



G^^\p,v) 



tE 



,n7iu. 



cos 



9 1 
9 1 



d'^-^p 



o^P-P 



{P + {^f + m^y 



L {2tiY- 



tE 

n=l 



cos(^^) / d p 



L 



oW-P 



For simplicity let us choose m = and p = 0. Using the fact that d p = p dpdVtd-i and 



{p + {^f + m^Y 



(42) 



(43) 



/ dVtfi-i = ^""d^i , a straightforward calculation yields 



.(2) 



Gq '{p,u,v)\p=o = Bi{d,L,u) + B2{d,L,v) + Bs{d,L,v) + B4{d,L,u) 



(44) 



where 



Bi{d,L,u) = — h2{d) / dkk"^ ^ coth Lk cosh ku. 



(45) 
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B2{d, L, v) = -j-hid) f dk k^~^ coth kL cosh kv. (46) 

Bsid, L, v) = j-h2{d) I dk k'^-^ sinh kv. (47) 

and finally 

B4{d,L,u) = -j-h2{d) f dkk'^-^ sinh ku. (48) 

It is worth mentioning that the structure of the divergences of Eqs.(^S|-^Sp are the same as for the 
tadpoles, as we will see. In the next section we will analyze the renormalization program for the 
two-point functions in both cases of DD and A^A^ boundary conditions. 



4 Analysis of the ultraviolet divergences oi Tjjjj[L^m^d^z 
and Tnn{Lj m, d, z 



The aim of this section is to analyze the structure of the divergences of the bare two-point functions 
for both cases DD and A^A^ boundary conditions. Let us start from the expression of the vacuum 
activity for the case of DD boundary conditions, i.e.. 



Using trigonometric identities and also the relation |24 



^ cosnx 1 7r cosha(7r — x) , , 



^^ 'n? + a? 2a? 2a sinh 



na 



which is vahd for < x < 2tt, it is easy to show that the vacuum activity in the case of DD b.c. 
is given by 

g 1 °° /■ 1 

TML, -, d, z) = ^j^^^.j:^ j d'-'P ^p^^n_.y^^2^ - 9fM L. -, z) (51) 

where 

f(T ^ \ 1 1 [ ^d-i 1 cosh((L-22:)(p^ + m^)^) 

f2(L,m,d,z) = -- — :-r-r / " P 1 -^ • (52) 

^ ^ 2{2nY-^J ^(p^ + ^2)| sinh(L(p2 + ^2)i) I ) 

In an analogous way, it is also possible to calculate the vacuum activity for the A^A^ b.c. i.e. 

T]\rj^{L, m, d, z) and we obtain 

o 1 °° /■ 1 

TNN{hm,d,z) = —^—^J2^J rf'^-ip______ + ^/2(L,m,d,z). (53) 

Since T£,j:,{L, m, d, z) and Tnn{L, m, d, z) have the same functional form, both have the same kind 
of ultraviolet divergences. Let us define fi{L,m,d) by: 

11°°/" 1 

The equation above has ultraviolet divergences, but it is (formally) proportional to the tadpole 
in finite temperature field theory, after the identification: /3 = 2L. To deal with the divergences 
of the one-loop two-point function at finite temperature we have to do frequency sums and [d — 
1) dimensional integrals for the continuum momenta. The most popular method to deal with 
Matsubara sums is to analytic extension away from the discrete complex energies down to real axis 



with the replacement of the energy sums by contour integrals p5[ . We prefer to use dimensional 



regularization in the continuum [^, and afterwards to analytically extend the modified Epstein 
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zeta function which appear after the dimensional regularization [^. Since the forniahsm has 
already been developed by Malbouisson and Svaiter in [^, we will only sketch the procedure 
here. First we have to use a well known result of dimensional regularization, i.e. 



Ui ^ 



and let us define the modified Epstein zeta function C(z, a) by: 



d'^k 7r2 d 1 

r(s- 0)3137, (55) 



00 -| 

az.a)= Y: . . ^ ^,y a'>0, (56) 



Jl=— 00 



which is analytic for Re{z) > ^ It is possible to analytic extend the modified Epstein zeta function 



where the integral representation is valid for Re{z) < 1, p8 



^ (^2_^^2)-._.l-2. 



a 



T{z -\) . . [^ {t^ - lydt 



V TT — ^, - + 4 sm nz I 
T[z) Ji 



(57) 



p2nat 1 

Using Eqs.(p5D and ( p?]) in Eq.(|^), we get a polar part (size independent) plus a size dependent 
analytic correction. It is clear that the mass counterterm generated by fi{L, m, d) is size indepen- 
dent, as the finite temperature field theory has no temperature dependent counterterm. The first 
interesting result of the paper is given by Eqs.(^) and and ([5^). The tadpole graphs expressed by 
Tdd{L, m, d, z) and Tnn{L, m, d, z) have the same z dependent part in modulus but with opposite 
signs. From the above discussion it is possible to understand the finiteness of the renormalized 
stress-tensor of an electromagnetic field near a flat prefectly conducting plate. Although the ex- 
pectation value of the squared electric and magnetic field are divergent, a delicate cancellation 
makes the renormalized stress-tensor finite. As the size-dependent parts of ToD{L,m,d^z) and 
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Tnn{L, m, d, z) have the same functional form and opposite signs, and recalhng that it is possible 
to treat the electric and magnetic modes separately (where one obeys DD and the other A^A^ b.c, 
we automatically obtain a finite result for the vacuum expectation value of the stress-tensor of the 
electromagnetic field. It is important to stress that when the conducting boundary is curved, the 
energy density diverges on the boundary [^ . 



To shall deal with the renormalization program in the one-loop approximation in the next 
section, also discussing, for the sake of completeness, the issue of IR divergencies in different 
numbers of dimensions. 

5 Analysis of the ultraviolet and infrared divergences of 
the 2;-dependent part of the tadpoles 

We will again use the fact that (i'^~^p = p'^~^dpdQd-i and / dfld-i = It^tt- It should be noted 

that, had we chosen m^ = 0, the ultraviolet divergences would have kept the same polar structure. 

Consequently, for simplicity let us choose again m = 0, and for reasons that will become evident 

latter, we first assume d > 3. The special case rf = 3 is discussed at the end of this section. 

Defining h2{d) by: 

111 

it is possible to write f2{L,m,d,z)\m=o as 

1 f°° 

f2(L,m,d, z)\m=o = -h2{d) / dkk'^~^ coth. kL cosh 2kz 

2 Jo 
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h2(d) / dhk'^"^ cosh kzsinhkz. (59) 

JO 



In a general way, the regularization process is achieved introducing exponential cut-off regulators 
and after this the identification of the poles of the regularized quantity by means of the Laurent 
series expansion around some point i.e. the negative power portion of such series. Note that 
instead of imposing renormalization conditions over the IPI correlation functions we can simply 
subtract the singular part of the Laurent series around some point, by the introduction of the 
counterterms. Let us assume 2; 7^ and z y^ L. A straightforward calculation gives 



f2{L,m,d,z)\m=(} = -h2{d) 



/ dA;fc'^~^(cothA;L- l)cosh2fc^ 
Jo 



POO 

+ dk k'^'^{cosh 2kz - sinh 2kz) 

Jo 



(60) 



In the first integral for large z, {cothkL — 1) has the behavior: {cothkL — 1) ~ e""^^^. Moreover, 
the second integral in the above equation is ultraviolet finite for z ^ 0. Let us define x = kL and 
q = kz in the first and second integrals above respectively. Then Eq.(|60D becomes: 

1 1 r°° ■ 2z 

f2{L,m,d,z)\m=o = 2^^^^'>Td^ J c?a;x°'"^(cotha; - 1) cosh(— x) 

1 1 /""^ 

+ o^2(rf)-r^ / dqq'^-^{cosh2q-smh2q). (61) 

2 z"- ^ Jo 

The second term in the above equation gives us the well known result that for a massless minimal 
coupled scalar field < ip'^{x) > diverges as ^ if we approach the plate |^. In order to analyze 
the polar part of f2{L, 0, d, z), we use the definition of the Gamma function. Let us define /i(z/, /i) 



22 



and J2(/i, /?) by 



/•oo T 

/i(/i,z/)=/ dxx^'-^e-"'' = —T{fi), Re{fi)>0, Re{u)>0 
Jo z/^ 



(62) 



and 



/5 



J2(/i,/3)=/ t/a;a;''-^e-^^(cothx-l) = 2i-^r(^)C(/i,- + l) i?e(/3) > 0, i?e(/i) > 1, (63) 
JO 2 



where (^(z, a) is the Riemman zeta function defined by 

1 



C{z,a) = J2 



n=0 



n + a) 



Re{z)>l, a ^0,-1, -2... 



(64) 



Then, using Eqs. (^2]), ( |63D and (|6^) in Eq. (l6l|) we have that 



f2{L,m,d,z)\m=o 



^^^(ci)-^ 



2^~'^r(ci - 2) C(rf - 2, - + 1) + C(rf - 2, -- + 1) 



+ 



-^2(c^)r(rf-2). 



(2^)'^-2 
Using the definition of the zeta function, it is evident that: 



(65) 



1 ^ 1 1 1 ^ 

2^ / xc(-2 + /r ^Nd-2 ^" rd-2 2^ 



d-2- 



(66) 



^'-\=o(n + (l + f ))''-' (L-.)^-2 ^'^-^.^i(n+(l-f); 
We see that the regularized f2{L, 0, d, z) has two poles of order (c? — 2) in z = and m z = L. 
Note that the residues of the poles in ^ = and in 2; = L are L-independent. Since the domain 
of analyticity of the zeta function is d > 3, the case d = 3 must be studied separately. Different 
treatments for d = ?> and d = A simply express the fact that infrared divergences are more severe 
in lower dimensions. 
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We will go back to Eg. (|52D , studying the case m^ 7^ 0, to see how the IR divergences pop 
up in the m^ —>■ limit. It is important to stress that, only in the N-N B.C.'s case we have IR 



divergences for massless fields, coming from the term n = 0, i.e., equations (|5lD and (|5^ are IR 
finite for m = 0. A straightforward calculation yields 



r rr J \ ^t.^j\f°°j P'^^'^ cosh((L-2z)(p2 + m2)5' 
f2{L,m,d,z) = -h2{d) / dp-^ -T )\ ^'^ ,/ • 

Defining cr = (p^ + 111^)2 and using the fact that d = 3, we have: 



(67) 



f2{L,m,d,z)\d=3 



M^) 



POO 

/ (icT (coth crL — 1) cosh 2(72; 

Jm 



(■00 

+ da (cosh 2az — sinh 2az) 

Jm 



(68) 



The second integral in the above expression is convergent for 2; 7^ 0, and defining v = 2a z, it 

becomes: 

1 r°° 1 

— /i2(3) / dve~" = —h2{3)T{l,2mz), (69) 

4:Z J2mz Az 

where r(a, x) is the incomplete gamma function. Consequently, we have a simple pole for z = 0. 
Again, the residue of this pole is L-independent. To complete the regularization procedure we 
have now to analyze the first integral of Eq. (|68|) : 



-/i2(3) / rfa(coth(TL-l)cosh2az = — /i2(3) / du + — /i2(3) / du . (70) 

2 Jm AL JmL e" — 1 AL JmL C" — 1 

The second integral in the right side of Eq. (^) is convergent and the first one has a simple pole 
at z = L, again with an L- independent residue. 



24 



From the discussion above, we can conclude that, in order to ehminate the uhraviolet diver- 
gences of the theory we have to introduce counterterms as surface interactions, and consequently 
the full action will have the following form for both fields (fi and ip2'. 

S{^) = j\z j d^'hi^id^f + ^mV' + ^A^^) + I d'-hic,^\f,0) + c,^'{f,L)). (71) 

As we have already taken care of the ultraviolet divergences, let us study the infrared divergent 
piece (for m = 0) of f2{L,m,d,z)\d=3- Let us call this piece f2{L,m,d,z)\d=3- Note that we 
introduce an ultraviolet cut-off in order to use the Bernoulli representation of the integrand. 

f;iL,m,d,z)\d=3 = jYh2iS)( du^^—+ du^^-] . (72) 

'±L/ \JmL e — i JmL C — i / 

Writing the integrand using the Bernoulli polynomials it is not difficult to show that 

1 z 27r 

f2{L,m,3,z) = — /i2(3)5o(-)/n(— ) + regular part {f;{L,m,z)). (73) 

When m ^ we have a logarithmic divergence which is z dependent. 

Going back to the case of the ultraviolet divergence, some authors claimed that the introduction 
of surface counterterms is against the spirit of the renormalization program. In our case, however, 
it is possible to change the model, by adding a new interaction term, in such a way that the 
ultraviolet divergences coming from the f2{L, m, d, z) contributions corresponding to each field are 
compensated. One possibility is to consider the 0(2) symmetric model with a ^\'^\ interaction 
term. Then the z-dependent part of each tadpole cancel each other out. The situation is similar to 
the case of supersymmetric theories, where the finiteness of some correlation functions is achieved 
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by a balance between bosonic and ferniionic loops. The 0(2) symmetric model with an interaction 
term ^pl^p"^ {DD and A^A^ B.C.'s) developes a size dependent mass Am^ proportional to gL'"^, as 
the (f"^ model at finite temperature. In the same way as temperature can solve the IR problem in 
some QFT models, finite size effects can also cure these divergences, when a resummation can be 
implemented. 

As stressed by many authors, in the electromagnetic case, the origin of the unboundedness 
renormalized stress-tensor near a curved surface has the origin in the unphysical nature of classical 
"perfect conductor" boundary condition. Let us suppose the following physical situation. For the 
low energy modes the manifold is [0, L] x ^'^^^ and for the high energy modes we have S^ x 3?^^^^, 
i.e. let us assume a sharp cut-off and for fc„ < A we have DD boundary conditions and for 
/cn > A we have periodic B.C.'s. The high frequencies do not "see" the mirror a.t z = 0, and 
translational invariance is maintained only for these modes. If the collapse of the renormalization 
program (removing infinities from perturbative calculations using only "bulk" counterterms) is 
related with the break of translational invariance, our improved model must be renormalizable. A 
further study of this model may be of interest. A different possibility is to construct an effective 



action for the slow-modes and after this imposing the DD or A^A^ b.c. [pT|||32 | 



6 Conclusions 

In this paper we studied finite size effects in an interacting field theory, with broken translation 
invariance. We calculated the vacuum activity for an anisotropic model, between two parallel plates 
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in ad dimensional Euclidean space. It has been possible to obtain closed expressions for < fHx) > 
and < v^2(^) >5 for fields satisfying DD and NN B.C.'s, respectively. For different shapes, 
the complexity and the number of technical difficulties increase enormously, but the multiple 
scattering method can be used in these cases [0], at least for small curvatures of the boundaries. 
We presented a model having the interesting property that the 2;-dependent part of the tadpole 
graphs for DD and A^A^ B.C.'s have the same modulus and opposite signs. This fact could explain 
the boundedness of the renormalized vacuum expectation value of the energy-stress tensor of the 
electromagnetic field in the Casimir-like configuration. 

We have also seen that, to renormalize the theory, counterterms corresponding to surface 
interactions are required. One can, however, avoid this difficulty by equipping the model with a 
(flfl interaction. Then, the z dependent pieces of each tadpole cancel each other out, and the 
two fields develop a size dependent mass Am"^ proportional to gL~'^, as for the single cp^ model at 
finite temperature. 

There are several directions in which the finite size effects for systems with breaking of trans- 
lational invariance which may deserve further research. To mention a few of them: the study 
of interacting fermions, the non-linear a model in domains with one finite direction and {d — 1) 



infinite directions ||3^, and finally as a straightforward extension of this work, the study of the 
0(2) symmetric model at the two- loop approximation. 

As discussed in the Introduction, one should prove that the the renormalization program can 
be implemented beyond the one-loop approximation, where overlapping divergences emerge. 
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